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It is known that the free surface of an axisymmetric viscous film flowing down the outside of a 
thin vertical fiber under the influence of gravity becomes unstable to interfacial perturbations. We 
present an experimental study using fluids with different densities, surface tensions and viscosities to 
investigate the growth and dynamics of these interfacial perturbations and to test the assumptions 
made by previous authors. We find the initial perturbation growth is exponential followed by a slower 
phase as the amplitude and wavelength saturate in size. Measurements of the perturbation growth 
for experiments conducted at low and moderate Reynolds numbers are compared to theoretical 
predictions developed from linear stability theory. Excellent agreement is found between predictions 
from a long-wave Stokes flow model (Craster & Matar, J. Fluid Mech. 553, 85 (2006)) and data, 
while fair agreement is found between predictions from a moderate Reynolds number model (Sisoev 
et al, Chem. Eng. Sci. 61, 7279 (2006)) and data. Furthermore, we find that a known transition in 
the longer-time perturbation dynamics from unsteady to steady behavior at a critical flow rate, Q c , 
is correlated to a transition in the rate at which perturbations naturally form along the fiber. For 
Q < Qc (steady case), the rate of perturbation formation is constant. As a result the position along 
the fiber where perturbations form is nearly fixed, and the spacing between consecutive perturbations 
remains constant as they travel 2 m down the fiber. For Q > Q c (unsteady case), the rate of 
perturbation formation is modulated. As a result the position along the fiber where perturbations 
form oscillates irregularly, and the initial speed and spacing between perturbations varies resulting 
in the coalescence of neighboring perturbations further down the fiber. 

PACS numbers: 47.20.-k, 47.20.Dr, 47.55.df, 47.85.mb 



I. INTRODUCTION 

Coatings are commonly applied to the exterior of thin 
cylindrical wires or fibers to provide protection and/or 
enhance performance (e.g. electrical wire and fiber-optic 
cable). Methods of coating include extruding a fiber 
through a die (die coating) or drawing a fiber from a 
liquid bath (dip coating) [E 0, & H H . During the coat- 
ing process, a uniform liquid film can become unstable 
to interfacial perturbations that may develop further into 
droplets [1,0]. This effect, which detracts from the qual- 
ity of a coating, has inspired a wide array of studies on 
the formation and motion of perturbations on cylindrical 
fibers [|, 0, 1 & M EH, El, El • 

Fibers can also be coated by a continuously-fed ax- 
isymmetric fluid flow down the length of a vertical fiber 
(see Figure [T|) as has been examined in several analytical 
and experimental studies [3, [H [M E3, H El- The 
geometry of the unperturbed flow is an annular film with 
a fixed internal boundary and a free surface at the outer 
fluid-air interface. It is well known that the free sur- 
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face of this annular film becomes unstable to interfacial 
perturbations, as shown in Figure [T] Herein we present 
an experimental study on an annular viscous film with a 
particular focus on the initial formation and longer-time 
dynamics of perturbations along the film free surface. 




FIG. 1: An annular film of viscous fluid flowing down the 
outside of a thin nylon fiber; the surrounding fluid is air. 
Note that the figure has been rotated by 90° with gravity 
acting towards the right. Perturbations develop along the 
free surface some distance down the fiber; once formed, these 
perturbations continue to travel down the fiber. Image length 
= 9.7 cm. 

A related problem to annular films is on the motion 
of cylindrical jets, which in contrast have no fixed in- 
ternal boundary. Analytical studies on the motion and 
stability of inviscid and viscous jets dates back to the 
work of Plateau [H, Lord Rayleigh [HH1], Weber [H 
and Chandrasekhar [24j]. It is known that capillary ef- 
fects drive perturbation growth along the jet free sur- 
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face, often referred to as the Plateau-Rayleigh instability. 
Analytical results developed from temporal linear stabil- 
ity theory 2j| w ere tested in experiments by Donnelly 
& Glaberson 25|, who found strong (fair) quantitative 
agreement between the theoretical and measured disper- 
sion relation for an inviscid (viscous) jet. The dynamics 
of free surface perturbations along cylindrical jets and 
annular films are, however, quite different. In the cylin- 
drical case jet breakup occurs when the perturbations 
become sufficiently large (25[, while in the annular case 
the large amplitude perturbations remain connected by 
a liquid film [11 [H 111]. 

Recently, several theoretical studies have analyzed the 
temporal linear stability of an annular viscous film flow- 
ing down a vertical fiber in the Stokes [3, EB] and mod- 
erate Reynolds number [l?], EH limits; the base flow in 
these studies is assumed to be a steady, unidirectional 
parallel flow. Here we test these results by determining 
whether: (i) the base flow used in [3, EE EE EH matches 
the experimental flow; and (ii) the dispersion relations 
derived in the Stokes and moderate Reynolds number 
limits [HI, EH correctly predict the nascent growth of 
perturbations measured in low and moderate Reynolds 
number flows. 

As free surface perturbations travel down a vertical 
fiber, many interesting phenomena occur [ljL EM El> Ell • 
In experiments, Kliakhandler, Davis and Bankoff (KDB) 
observed three types of behavior far down the length of 
the fiber (« 2 m) [l4j]. At the highest flow rate (regime 
a), the film between perturbations is thick and uniform, 
and faster moving perturbations collide into slower mov- 
ing perturbations (unsteady behavior) . At an intermedi- 
ate flow rate (regime b), the spacing, size and speed of 
the perturbations is constant so that no collisions occur 
(steady behavior). And at the lowest flow rate (regime 
c), the fluid periodically drips from the tank, rather than 
jets as with the higher flow rates, creating a regular spac- 
ing between perturbations near the tank outlet. The long 
time between drips allows the film connecting consecutive 
perturbations to thin and subsequently become unstable 
to smaller capillary perturbations. (Figure 1 in Jl4j il- 
lustrates these three regimes of behavior.) Simulations 
of a Stokes flow model developed by KDB qualitatively 
captured two of the three observed behaviors (regimes b 
and c), while the behavior associated with the highest 
flow rate (regime a) could not be replicated [T^J. 

Craster & Matar (CM) [l5| select a different scaling 
than KDB to derive an evolution equation for the free 
surface. Using traveling wave solutions, their Stokes flow 
model quantitatively predicted the perturbation speed 
and height of regime a measured by KDB. The model 
also qualitatively captured regime c, though the steady 
pattern of perturbation spacing found in regime b could 
not be matched with traveling wave solutions. In experi- 
ments, CM observed regime b near the tank outlet, how- 
ever, they found the regularly spaced pattern of pertur- 
bations disassembled itself further down the fiber. From 
this observation, CM concluded regime b is a transient 



rather than steady regime [15j |. 

The contradiction in observations of regime b (steady 
behavior) by KDB and CM motivated us to look more 
closely at the steady and unsteady states by examining 
the dynamics of the perturbations where they initially 
form along the fiber. In experiments using fluids with 
different densities, surface tensions and viscosities we ob- 
serve regimes a (unsteady), b (steady) and c (dripping) 
as described by KDB [l4[ , though the focus of this paper 
is on regimes a and b. In our experiments, we observe 
the flow transitions abruptly from unsteady to steady 
behavior at a critical flow rate Qc (the value of Q c is 
dependent on the particular fl uid) [26| . In a recent inde- 
pendent study, Duprat et al. [1 91 ] explain the transition 
from regime a (unsteady behavior) to regime b (steady 
behavior) as a transition from a convective to absolute 
instability. In their experiments with silicone oil using a 
range of fiber and orifice radii, they find the transition 
occurs only at intermediate film thicknesses and for suf- 
ficiently small fiber radii; at thin or thick film thickness, 
they find the perturbation behavior remains convective 
(unsteady) [19| . These criteria may explain why CM did 
not observe the steady dynamics of regime b in their ex- 
periments with silicone oil. Here we find the transition 
from unsteady to steady behavior is also correlated to 
the rate at which perturbations naturally form along the 
fiber. For Q < Q c (steady case), the rate of pertur- 
bation formation is constant. As a result the position 
along the fiber where perturbations form is nearly fixed, 
and the spacing between consecutive perturbations re- 
mains constant as they travel 2 m down the fiber. For 
Q > Qc (unsteady case) , the rate of perturbation forma- 
tion is modulated. As a result the position along the fiber 
where perturbations form oscillates irregularly, and the 
initial speed and spacing between perturbations varies 
resulting in the coalescence of neighboring perturbations 
further down the fiber. 

The paper is organized as follows. The experimental 
setup and properties of the unperturbed flow are pre- 
sented in Section [Til Measurements of the perturbation 
growth are compared to analytical predictions for Stokes 
and moderate Reynolds number conditions in Sect ion HTT1 
The perturbation behavior exhibited in regimes a (un- 
steady) and b (steady) near the tank outlet is closely 
examined in Section HVl Conclusions are provided in 
Section N\ 



II. EXPERIMENTAL SETUP & DETAILS 



A. Experimental Apparatus 

The experimental setup consists of viscous fluids, a 
reservoir tank/orifice assembly, nylon fishing line, a high- 
speed digital imaging camera, illumination, a computer 
and edge-detection software (a schematic of the experi- 
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(a) 



Reservoir Tank 



Nylon Fiber 



(b) 



C.mer. □ 




TABLE I: Fluid properties and experimental conditions. 



, r-S(z.t) 



FIG. 2: (a) Schematic of experimental setup, (b) Cross- 
section of an annular film flowing down the outside of a fiber 
of radius r/ (not to scale). The unperturbed film radius 
measured from the fiber centerline is ro, unperturbed film 
thickness ho, perturbation amplitude T)(z,t) and perturbation 
wavelength A. 



ment is shown in Figure a)). The reservoir tank (6 L 
capacity) is graduated at 100 mL increments to measure 
flow rate. An orifice, machined with a flat edge (inner ra- 
dius = 0.11 cm, outer radius = 0.16 cm, length = 2.6 cm), 
is attached to the bottom of the tank to ensure a repro- 
ducible solid/fluid/air contact line in the experiment. A 
nylon fiber (rt = 0.029 cm) anchored from above, passes 
through the center of the tank/orifice assembly, and is 
held vertically plumb with weights attached 2 m below 
the orifice. The fluid, which is gravitationally forced from 
the tank, coats the fiber to create an annular film. To re- 
duce air currents and other noise during data collection, 
the entire apparatus is enclosed by an aluminum frame 
with plastic sheet sidewalls and top. 

The motion of the annular film is recorded using a high- 
speed digital imaging camera (Phantom v4.2) at rates be- 
tween 1000-4000 frames/s and an image size of 64 x 512 
pixels 2 with the camera focused on approximately the 
upper 10 cm of the fiber. Illumination is obtained us- 
ing silhouette photography following [27], with a 250 W 
lamp, an experimental grade one-way transparent mirror 
(Edmund Scientific, A40,047) and high contrast reflec- 
tive screen material (Scotchlite 3M 7615). Movies of the 
annular film are recorded and downloaded to a computer 
using camera software. The free surface of the film is 
determined from movie images using an edge-detection 
algorithm. The algorithm locates the free surface by in- 
terpolating the maxima positions in a gradient image; the 
gradient image is produced using the Frei-Chen operator 
|28| . The algorithm can detect the edge of the film free 
surface to within approximately l/10 th of a pixel which 
for the screen resolution in our experiments corresponds 
to w 0.002 cm. 

The experimental fluids consist of castor oil, vegetable 
oil (Crisco) and an 80:20 glycerol/ water solution (by 
weight). The temperature, density (p), surface tension 
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Castor Oil 


21.9 


0.94 


36.8 


8.48 


1000 


0.0190 


Vegetable Oil 


21.6 


0.92 


34.3 


0.58 


4000 


0.0180 


Glycerol Solution 


21.2 


1.21 


60.4 


0.54 


2000 


0.0186 



(7), dynamic viscosity (77) of the fluids, and the fram- 
ing rate and screen resolution used in the experiments 
are listed in Table |TJ The surface tension was measured 
at room temperature using a Fisher 21 tensiomat and 
viscosity was measured using a temperature controlled 
cone and plate rheometer (Brookfield, Model DV-III+). 
The fluid temperature varied by less than 4.6%, 1.9% 
and 1.4% in the castor oil, vegetable oil and glycerol so- 
lution experimental runs, respectively. We note that the 
selection of fluids allows us to independently probe the 
influence of surface tension or viscosity on flow behavior 
since castor oil and vegetable oil have comparable surface 
tension while vegetable oil and the glycerol solution have 
comparable viscosity. 



B. Base Flow Properties 

In each experiment, the reservoir tank drained under 
the influence of gravity and the elapsed time and tank 
volume were recorded as the fluid passed each 100 mL 
mark to determine the flow rate (Q). Data was col- 
lected only while the flow was jetting from the orifice 
and an unperturbed region of the film free surface was 
present near the orifice (corresponding to regimes a and 
b [14(). Measurements of the flow rate during each ex- 
perimental run for vegetable oil (square), glycerol solu- 
tion (triangle) and castor oil (circle), is shown in Fig- 
ure [3]Ja). (Note that in each experiment the flow rate 
decreases as the tank volume decreases.) Figure [Ha) 
shows that the flow rate increases linearly over the range 
of tank volume used in each experimental run. Further- 
more, the flow rate decreases with increasing viscosity 
(Qcastor is an order of magnitude less than Q g i yC eroi and 
Q vegetable), and increases with increasing surface tension 

(Q glycerol Q vegetable) ■ 

The unperturbed film radius ro, measured from the 
fiber centerline to the unperturbed free surface (see Fig- 
ure Hlh)), was measured using edge-detection software. 
Figure [3jb) shows the unperturbed radius increases with 
flow rate in each experimental run. This trend is simi- 
lar to behavior observed in dip coating in which the film 
thickness increases with the velocity at which the fiber is 
drawn from the fluid source at sufficiently low withdrawal 
rates [l|, [3, H, @ . Comparing data for the glycerol solu- 
tion (triangle) and vegetable oil (square) experiments, we 
find that at a fixed flow rate higher surface tension (glyc- 
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FIG. 3: (a) Experimental flow rate as a function of the reser- 
voir tank volume, (b) Flow rate as a function of the unper- 
turbed film radius: experimental data (symbols), theory given 
by ((4]) with r/ = 0.029 cm (dotted); bars represent the resolu- 
tion of experimental measurements. Symbols correspond to: 
vegetable oil (square), glycerol solution (triangle) and castor 
oil (circle). 



erol solution) results in a thinner unperturbed film. In 
our experiments, the range of the ratio of unperturbed 
film thickness to fiber radius is 2.5 < h^/rj < 3.3, which 
means the films are thick. 

The relevant forces that characterize the dynamics of 
an annular film can be determined by considering two 
dimensionless groups, the Reynolds number and Bond 
number. Following CM we define characteristic 

length and velocity scales as 



£ = l/(P9ro), V = pgrl/n, 



(1) 



where g is the gravitational constant of acceleration, the 
Reynolds number as Re = pVC/p, which compares in- 
ertial to viscous effects, and the Bond number as Bo = 
P9 r o 1 7j_which compares gravitational to surface tension 
effects [15]. The values of Re and Bo for our experiments, 
listed in Table HT1 indicate inertial effects are negligible for 
castor oil (Re ~ O(10 -2 )), in contrast to the vegetable 
oil and glycerol solution experiments (Re ~ O(10)), and 
surface tension dominates over gravitational effects in the 
formation of interfacial perturbations for all three fluids 
(Bo-OtlO- 1 )). 



TABLE II: The Reynolds number and Bond number measured 
in experiments. 



Fluid 


Re 


Bo 


Castor Oil 


0.049 - 0.054 


0.26 - 0.32 


Vegetable Oil 


9.4 - 11.6 


0.26 - 0.40 


Glycerol Solution 


27.2 - 31.2 


0.23 - 0.30 



Experimental [lj, [15J and theoretical [lj, |l5j, ll7|, |l8| 
studies often assume the base flow of an annular film 
is steady, unidirectional and parallel. Under these con- 
ditions, the unperturbed flow with free surface located 
at r = ro and constant pressure field p(r, z, t) — po is 
described by the boundary value problem for the axial 
velocity w(r) 

pr^ 1 d r (rd r w) + pq = 0, iv(rf) = 0, d r w(ro) — 0, 

(2) 

where the boundary conditions include no-slip at the 
fiber and zero tangential stress at the free surface. 
Eqs. ([2]) can be solved exactly for the axial velocity to 
obtain 



w(r) 



pg_ 

4,, 



2r 2 \n(r/r f )+r 2 f ~r 2 



(3) 



Using ([3]), the flow rate of an annular film can be ex- 
pressed in terms of rj and rp 



Q{rf,r ) 



2n 
8p 



rw(r) dr 

(Ary a +44\n(r /r f )-rj-3r 4 )(A) 



By comparing ([J) to experimental data, we can test the 
assumption that the base flow of an annular film is a 
steady, unidirectional parallel flow. 

Figure [21b) shows a comparison of the flow rate as a 
function of the unperturbed film radius for vegetable oil, 
glycerol solution and castor oil measured directly in ex- 
periments (symbols) and using ^ with rf = 0.029 cm 
(dotted line). We find excellent agreement between ^ 
and the experimental data in the castor oil (circle) and 
vegetable oil (square) experiments which indicates the 
flow in these experiments is well approximated by a unidi- 
rectional parallel flow. The theory overestimates Q by as 
much as 12% in the glycerol solution experiment (trian- 
gle). This is not surprising since the assumptions on the 
flow are equivalent to steady Stokes flow. With Re w 30, 
we suspect the Stokes flow assumption is not valid for 
the glycerol solution experiment. Based on the compar- 
ison in Figure [3][b), we estimate ((3]) and ([4]) are valid 
when Re < 10. The experiments conducted by KDB [3| 
and CM [l5[ meet this criterion, thus we conclude their 
assumption that the flow is steady, unidirectional and 
parallel is indeed valid. 

In the proceeding sections, we examine the dynamics 
of a perturbed annular film including the initial forma- 
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tion and longer-time dynamics of interfacial perturba- 
tions along the free surface. 



Reynolds number flow [l8j limits. 



III. PERTURBATION GROWTH 



A. Linear Stability Results 



The image in Figure [T] illustrates the capillary insta- 
bility an annular viscous film undergoes as the unper- 
turbed free surface becomes unstable to undulations that 
develop into large amplitude perturbations. Next we 
present experimental observations on the growth of these 
interfacial perturbations and compare their initial growth 
to theoretical predictions developed from linear stability 
analysis. Before proceeding, we first recount relevant sta- 
bility results developed in the Stokes [l5| and moderate 
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1. Stokes Flow 

Craster & Matar [l5[ derive a long-wave Stokes flow 
evolution equation for the free surface of the annular film, 
r = S(z,t), under the assumption that the unperturbed 
film radius ro is small relative to the capillary length 
l c = l/pgro (i.e., Bo = r /l c = pgr oh < 1) and the 
Reynolds number is sufficiently small (Re < 0(1)), to 
obtain 



n, ( j - (Bo) 2 dz,S ) I 



x (a 4 - AS 2 a 2 + 3S 4 - 45 4 ln(5/a)) 
I 



= 0, 



(5) 



where S, z, t are dimensionless variables satisfying the 
scalings 



S = r S, 



Cz, t = Ct/V, 



(6) 



and a = r//fo- Conducting a linear stability analysis by 
perturbing about the base flow 



S(z,t) = 1 + S 1 e^ m£+ai \ 



(7) 



where m is the (real) dimensionless wavenumber and a is 
the (complex) dimensionless growth rate CM obtain the 
following dispersion relation for the growth rate 



where m is the (real) dimensionless wavenumber and a 
is the (complex) dimensionless growth rate, Sisoev et al. 
obtain a dispersion relation for a satisfying 

<j 2 + (ia 1 ^-a 1 ^m)<j+- (— a ,3»™ 3 — a,o,im+ia Q . ) = 0, 

+ £ . (10) 
where e = ho/rf and the constant coefficients 

ai,o, a i,ii a o,3i 0-0,1, a a,o arc defined in the Appendix. The 

variables h, z, t are dimensionless quantities satisfying the 

scalings 
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(BoW - l)(a 4 



4a 1 



3 + 41n(a)) 



im 



(a 2 - 1 - 2 ln(a)) 



(8) 



hnh. 



where 



kU 



(11) 



nolds Number Flow 



2. Moderate 



In an analytical study, Trifonov [l7l | derived model 
equations for fluid flowing down the inside or outside 
of a vertical cylinder at moderate Reynolds number; the 
model includes evolution equations for the film thickness, 
h(z,t) — S(z,t) ~ Tf, and volumetric flow rate, q(z,t). 
In a recent study, Sisoev at al. [II] rescale Trifonov's 
equations for flow down the outside of a vertical cylin- 
der, casting the model in terms of a gen eralized falling 
film model [H] (see Eqs. (11)-(13) in [18]). Conducting 
a linear stability analysis of the rescaled equations by 
perturbing about the base solution 



h(z,t) = l + hie*™*-**), 
q(z,t) = l + fce'f"**-^, 



(9a) 
(9b) 



P9K 

7 



1/3 



u 



2nrfho 



(12) 



represent a stretching parameter and characteristic veloc- 
ity scale, respectively. We note that the long-wave model 
used by Sisoev et al. was derived under the assumption 
that k 2 < 1. 

The stability results developed by Craster & Matar 
[III and Sisoev at al. [l8[ are temporal analysis (since 
to 6 K and <r e C), and thus model the case in which 
interfacial perturbations grow in amplitude everywhere 
along the film [3(| • Our interest is in testing these stabil- 
ity predictions by comparing the theoretical dispersion 
relations ([5]) and (jTUJ) to the growth rate of perturba- 
tions measured in experiments conducted in the Stokes 
and moderate Reynolds number flow limits. 
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B. Experimental Observations of Perturbation 
Formation 




(a) (b) (c) (d) (e) (f) (g) (h) (i) (j") 

FIG. 4: An annular film of castor oil flowing down the outside 
of a thin nylon fiber. The film loses uniformity approximately 
5.4 cm from the orifice indicated by the circled area in frame 
(b). Subsequent frames track the position and growth of this 
nascent perturbation. Time between images is 0.42 s and 
elapsed time = 3.7 s. Image height = 9.7 cm. 



Figure 2] shows a series of images tracking the forma- 
tion of a perturbation along an annular film of castor 
oil. In frame (b), a small amplitude perturbation first 
appears along the film approximately 5.4 cm from the 
orifice (circled region). Frames (c)-(j) track the position 
of this perturbation as it grows in amplitude and satu- 
rates in size. Once formed, the perturbation continues 
moving down the fiber (not shown). Since the pertur- 
bation grows in amplitude as it travels down the fiber, 
the flow is spatially (convectively) unstable rather than 
temporally (absolutely) unstable to perturbations [3(| ■ 

To characterize the growth of a perturbation we mea- 
sure the amplitude ?y (half the radial distance from first 
minima to first maxima) and the wavelength A (the ax- 
ial distance from first to second maxima) as shown in 
Figure Hfb) using edge-detection software; both measure- 
ments are made in the moving reference frame of the per- 
turbation. The data shown in Figure[5]corresponds to the 
perturbation followed in Figured] Figure Ufa) shows the 
nascent growth of the amplitude is exponential (inset) 
followed by a slower phase as the perturbation saturates 
in size (fjsaturate = 0.035 cm). The growth rate for the 
initial formation of the amplitude is determined from a 
least squares fit of the data to an exponential function 
yielding the dimensional growth rate Udim — 0.959 s _1 
(fit indicated by dotted line in Figure Ufa)). The wave- 
length of this perturbation decreases from A = 0.98 cm 
to 0.80 cm during the time interval that the amplitude 
grows exponentially (0 < t < 2.34 s), before saturating 
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0.05 r 



0.025 




(a) 




FIG. 5: (a) Amplitude (77) and (b) wavelength (A) of a 
nascent perturbation as a function of time for castor oil 
(Q — 0.0194 cm 3 /s); experimental data corresponds to the 
perturbation tracked in Figure [4] (a) The initial amplitude 
growth of the perturbation is exponential (shown in inset on 
log-linear scale) followed by nonlinear saturation. Dotted line: 
fit of data to r)(t) = 0.0023e°' 959t corresponding to growth 
between frames (b)-(g) in Figure [3] (b) The wavelength de- 
creases during the time interval that the amplitude grows ex- 
ponentially and then saturates in length as the amplitude 
saturates in size. 



in length to A = 0.83 cm, as shown in Figure [5jb). The 
decrease in A during the exponential phase of growth indi- 
cates the annular film is unstable to a range of wavcnum- 
bers (= 27t/ A), rather than to one fixed value. 

The behavior displayed in Figure O for the amplitude 
and wavelength is typical of observations made in the cas- 
tor oil, vegetable oil and glycerol solution experiments. 
Comparison of perturbation growth to the Stokes © or 
moderate Reynolds number (|10[) dispersion relations de- 
pend on the flow conditions in each experiment, specif- 
ically on the Reynolds and Bond numbers (provided in 
Table HI)) . In the castor oil experiments, Bo ~ O(10 _1 ) 
and Re ~ O(10 -2 ), thus satisfying the requirements of 
the Stokes model (Bo « 1, Re < 0(1)). Since Re w 30 in 
the glycerol solution experiments, inertial effects cannot 
be ignored, and so we compare this case to the moder- 
ate Reynolds number model. With Bo ~ O(10 _1 ) and 
9.4 < Re < 11.6, the vegetable oil experiments are on 
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the border of the requirements for the Stokes model. In 
this case we compare the experimental data to both the 
Stokes and moderate Reynolds number dispersion rela- 
tions. 
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FIG. 6: (Color online) Dimensionless growth rate as a func- 
tion of dimensionless wavenumber for experiments with (a) 
castor oil and (b) vegetable oil. Circles correspond to ex- 
perimental data of the average growth rate and wavenumber 
measured over several perturbations at (a): Q = 0.0194 - or- 
ange, 0.0171 - black, 0.0148 - cyan, 0.0123 - red, 0.0111 - green 
and 0.0103 - blue (cm 3 /s); (b): Q = 0.416 - orange, 0.383 - 
black, 0.342 - cyan, 0.252 - red, 0.211 - green and 0.170 - blue 
(cm 3 /s). Vertical bars represent the standard deviation of a 
and horizontal bars represent the range of m measured dur- 
ing the period of exponential growth over all perturbations 
measured. Corresponding colored curves represent theoreti- 
cal prediction given by the Stokes flow dispersion relation (JSj 
of IZe(a) plotted over the range of m measured in experiments 
0. 



Figure [5] shows a comparison of the measured ampli- 
tude growth rate to the dispersion relation developed by 
Craster & Matar in the Stokes flow limit ((8]) for (a) cas- 
tor oil and (b) vegetable oil at various flow rates. At 
a given flow rate, the growth rate for several perturba- 
tions was measured (8-12 perturbations for castor oil, 
and 15-44 for vegetable oil). The average dimensionless 
growth rate (a = adi m C/V) and dimensionless wavenum- 
ber (m = 2t:C/X) measured over all the perturbations is 



denoted by a circle with each color corresponding to a dif- 
ferent flow rate. The vertical bars represent the standard 
deviation of all the growth rates measured at a given flow 
rate during the period of exponential growth. Since the 
perturbation wavelength decreases over a range of val- 
ues during the exponential phase of growth, we cannot 
assign a single wavenumber to its growth. Instead, the 
horizontal bars represent the range of wavenumber mea- 
sured during the period of exponential growth of all the 
perturbations. The corresponding colored curves repre- 
sent the real part of the growth rate predicted by (|8|) 
plotted over the range of wavenumber measured at each 
flow rate. We consider the theory to be in agreement 
with the experimental data (at a given flow rate) if the 
theoretical curve overlaps the rectangular region defined 
by the resolution bars of the data. Figure [5] shows that 
the Stokes theory is in agreement with four of the six cas- 
tor oil experiments and with five of the six vegetable oil 
experiments. In the other three experiments, the theory 
overestimates the measured values by 10 to 13%. The 
quantitative agreement between theory and experimen- 
tal data is excellent, a significant result considering: (i) 
the comparison is between a temporal stability theory 
and a spatial instability of the film, and (ii) the value 
of the Reynolds number in the vegetable oil experiments, 
Re w 10, is slightly higher than the criteria for the Stokes 
theory, Re < O(l). 

Figure [7] shows a comparison of the measured ampli- 
tude growth rate to the dispersion relation developed by 
Sisoev et al. in the moderate Reynolds number limit 
(fTU| for (a) vegetable oil and (b) glycerol solution at var- 
ious flow rates. The data and theory are presented in a 
similar fashion to Figure [5] with the exception that the 
dimensionless growth rate and wavenumber are given by 
g = o dimho I (kU) and m = 27t/io/(kA), and the growth 
rates for the glycerol solution experiments are averaged 
over 88 to 102 perturbations. We recall that the moder- 
ate Reynolds number model is valid as long as k 2 <C 1; in 
all of the experiments shown in Figure [71 k 2 ~ O(10 _1 ). 
We find the moderate Reynolds number model overesti- 
mates the measured growth rates by 28 to 50% in the 
vegetable oil experiments and 15 to 48 % in the glycerol 
solution experiments, as shown in Figure [7] Clearly, the 
Stokes model is more accurate at predicting the growth 
rate of the perturbations in the vegetable oil experiments 
than the moderate Reynolds number model. This is 
somewhat surprising since the vegetable oil experiments 
slightly exceed the Reynolds number limit of the Stokes 
model, but satisfy the assumption on k 2 for the moderate 
Reynolds number model. While the theoretical growth 
rates in Figure [7] are on the same order of magnitude 
as the measured values, the quantitative match between 
theory and data is not strong. 

We note that the range of the measured amplitude 
growth rate in experiments (indicated by the vertical bars 
in Figures [5] & [7J varies by fluid and flow rate. For cas- 
tor oil (Figure [Sfa)), the range is fairly small which we 
attribute to the low Reynolds number (Re ~ O(10 -2 )) 
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FIG. 7: (Color online) Dimensionless growth rate as a func- 
tion of dimensionless wavenumber for experiments with (a) 
vegetable oil and (b) glycerol solution. Circles correspond to 
experimental data of the average growth rate and wavenum- 
ber measured over several perturbations at (a): Q = 0.416 - 
orange, 0.383 - black, 0.342 - cyan, 0.252 - red, 0.211 - green 
and 0.170 - blue (cm 3 /s); (b): Q = 0.538 - orange, 0.493 - 
black, 0.437 - red, 0.381 - green and 0.325 - blue (cm 3 /s). 
Vertical bars represent the standard deviation of a and hori- 
zontal bars represent the range of m measured during the pe- 
riod of exponential growth over all perturbations measured. 
Corresponding colored curves represent theoretical prediction 
given by the moderate Reynolds number flow dispersion rela- 
tion (|10|) of Tm(a) plotted over the range of m measured in 
experiments [l8l ]. 



behavior in the blue glycerol solution data set is steady. 
(The notion of unsteady and steady perturbation behav- 
ior will be explained in detail in Section HVl ) Therefore, 
we find the range of amplitude growth rate of the pertur- 
bations is correlated to both the Reynolds number of the 
flow and the longer-time dynamics of the perturbations. 

Next, we examine the dynamics of perturbations after 
their initial formation and explain a physical mechanism 
that controls a known transition in the flow from un- 
steady to steady perturbation behavior. 



IV. STEADY AND UNSTEADY 
PERTURBATION DYNAMICS 

The dynamics of interfacial perturbations along an an- 
nular film flowing down a vertical fiber can be broken 
down into three essential stages: (i) initial exponential 
growth of the perturbation amplitude accompanied by a 
decrease in wavelength; (ii) nonlinear saturation of the 
perturbation amplitude and wavelength; and (iii) longer- 
time behavior in which the perturbation wavelength may 
(unsteady - see Figure [T]) or may not (steady) vary along 
the fiber; this last stage has been noted in other exper- 
imental studies [3, EH, Eil • Here we explain a physical 
mechanism that controls this third stage of dynamics. 
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in the experiments. For the experiments with vegetable 
oil (Figures [6]Jb) and Eta)) and glycerol solution (Fig- 
ure E£b)) with Q > 0.345 cm 3 /s, the range of the mea- 
sured amplitude growth rates is large. Naively, one could 
attribute this to the higher Reynolds number in these 
experiments (10 < Re < 30). This is, however, not the 
complete picture. Notice the range is significantly smaller 
for the glycerol solution experiment at Q — 0.325 cm 3 /s 
(blue (rightmost) data set in Figure [7fb)) . The difference 
in this data set compared to the other glycerol solution 
and vegetable oil sets is in the behavior of the perturba- 
tions. The perturbation behavior in the sets with a large 
range of amplitude growth rates is unsteady, while the 



FIG. 8: Coalescence of two perturbations along an annular 
film of glycerol solution. Q = 0.359 cm 3 /s, time between im- 
ages is 0.0125 s, elapsed time = 0.1 s, image height = 9.7 cm. 

In experiments with all three fluids, we observe the 
perturbation motion abruptly transitions from unsteady 
(regime a) to steady (regime b) behavior at a critical flow 
rate, Q c (Q c = 0.0095 cm 3 /s for castor oil, 0.119 cm 3 /s 
for vegetable oil and 0.345 cm 3 /s for glycerol solution) 
(26| , similar to observations made by Duprat et al. in 
their experiments with silicone oil Following KDB 
(l4| . we define the flow to be steady if no perturbations 
coalesce as they travel down the full length of the fiber 




FIG. 9: Space-time plot illustrating the motion of perturbations along a vertical fiber during (a) unsteady and (b) steady 
behavior. Lighter (darker) regions correspond to thicker (thinner) regions of the fluid interface with light characteristic lines 
indicating the location of perturbations along the fiber as a function of time. Experimental fluid: glycerol solution (Q c 



0.345 cm 3 /s), (a) Q = 0.347 cm 3 /s, (b) Q 
image is 0.58 cm below the orifice. 



0.336 cm /s, elapsed time = 8.09 s, image height = 8.22 cm. The top of each 



(« 2 m), and unsteady otherwise while the flow is jetting 
from the orifice. An example of unsteady behavior in 
which two perturbations coalesce is shown in Figure [5] 
Note: we will not be examining the dripping state, regime 
c [HI], which occurs at a lower flow rate, Qdrip < Qc- We 
find the transition from unsteady (Q > Q c ) to steady 
(Qdrip < Q < Q c ) behavior is robust in the sense that 
once an experiment transitions to steady behavior (as 
the flow rate decreases) it does not revert back to the 
unsteady state. 

The space-time plots in Figure[9]illustrate (a) unsteady 
and (b) steady perturbation behavior for experiments 
with glycerol solution. Each plot is focused on 8.22 cm of 
the film, with the top of each plot located 0.58 cm below 
the orifice; the time span of each plot is 8.09 s. The plots 
are created by mapping the radius of the free surface of 
the film, r = S(z, t), to a gray level with lighter (darker) 
gray level corresponding to thicker (thinner) regions of 
the free surface. The light characteristic lines in the 
plots indicate the location of perturbations as they move 
down the fiber, and their slope represent the speed of 
the perturbations. Two features in the space-time plots 
distinguish the unsteady and steady perturbation behav- 
ior. First, the location along the fiber that perturbations 
form (which we refer to as the boundary) oscillates irreg- 
ularly in the unsteady case and appears nearly fixed in 
the steady case [2f|. Second, in the unsteady case per- 
turbations coalesce as faster moving perturbations collide 
into slower moving perturbations (indicated by intersect- 
ing characteristic lines), whereas in the steady case per- 
turbations do not coalesce as they travel with the same 
terminal speed down the fiber (indicated by parallel char- 




FIG. 10: Space-time plot illustrating the motion of pertur- 
bations along a vertical fiber during unsteady behavior. Ex- 
perimental fluid: glycerol solution; Q — 0.516 cm 3 /s, elapsed 
time = 8.09 s, image height = 8.22 cm. The top of the image 
is 0.58 cm below the orifice. 



acteristic lines) [3, ■ The longer-time motion of the 
perturbations appears to be correlated to the motion of 
the boundary. Notice in Figure [5Ja) that large spatial 
variations in the boundary modulate the perturbation 
speed (i.e., the slope of the characteristic lines) which re- 
sults in coalescence events later down the fiber. In the 
steady case, there is no spatial variation in the boundary, 
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and as a result, the perturbations remain equally spaced 
as they travel with constant terminal speed down the full 
length of the fiber (not shown in Figure 02b)). Our ob- 
servations of the steady case (regime b) are consistent 
with those of KDB |l4[. Given the robustness of the 
steady dynamics in all of our experiments, we conclude 
this is not a transient state as CM report [HI]. Finally, 
we note that when the flow is unsteady, the oscillation 
frequency of the boundary increases with increasing flow 
rate; for example, compare the boundary frequency in 
Figures [1(a) andCTOl 
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FIG. 11: (a) The distance from the orifice (D) where pertur- 
bations initially form along the fiber (i.e., the location of the 
boundary) as a function of time, corresponding to the data 
shown in Figure[9ja) for the unsteady case at Q — 0.347 cm 3 /s 
(•) and in Figure EJb) for the steady case at Q = 0.336 cm 3 /s 
(+)• (b) Average distance (D aV g) from the orifice that per- 
turbations form as a function of flow rate Q for experiments 
with glycerol solution. Vertical bars represent the standard 
deviation of D over all the perturbations measured at a given 
flow rate. The dotted vertical line denotes the transition flow 
rate separating steady and unsteady perturbation behavior. 



To characterize the motion of the boundary, we mea- 
sure the distance from the orifice D that each perturba- 
tion initially forms along the fiber (at a fixed flow rate) 
using edge-detection software; a perturbation is detected 
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FIG. 12: Power spectrum of the distance from the orifice D 
that perturbations form as a function of time. Plots corre- 
spond to the data shown in Figure lllf a) in: (a) steady be- 
havior (+) at Q = 0.336 cm 3 /s and (b) unsteady behavior 
(•) at Q — 0.347cm 3 /s for glycerol solution. In the steady 
case, the fundamental frequency /o = 14.45 Hz, is the first 
harmonic in the power spectrum. In the unsteady case, the 
bandwidth supporting the fundamental peak is larger than 
the steady case. 



when its amplitude (?y) initially exceeds 1 / 10 th of a pixel 
or » 0.002 cm. The data shown in Figure flTTa) corre- 
sponds to the space-time plots of the unsteady (•) and 
steady (+) experiments shown in Figure [9l Figure fTTT b) 
is a plot of the average distance from the orifice that per- 
turbations form (-Davg) as a function of flow rate in exper- 
iments with glycerol solution. The vertical bars represent 
the standard deviation of D over all the perturbations 
measured at a fixed flow rate and the dotted vertical line 
represents the transition flow rate, Q c . The distance that 
perturbations form from the orifice increases monotoni- 
cally with increasing flow rate. In the steady case, at a 
given flow rate the distance is nearly constant, whereas 
in the unsteady case, the range of distance that perturba- 
tions form increases with increasing flow rate. These re- 
sults are consistent with experimental observations made 
by Duprat et al. [l9l |. 
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FIG. 13: (a) Fundamental frequency as a function of flow 
rate measured from power spectra of D versus t data in ex- 
periments with glycerol solution. The dotted vertical line de- 
notes the transition flow rate separating steady and unsteady 
perturbation behavior, (b) The (normalized) cumulative in- 
tegrated power measured within the support of the first peak 
between fi ow = 6.18 Hz and fhigh = 19.77 Hz for the spec- 
tra shown in Figure [12T b). The Interquartile Region (IQR) 
is the frequency bandwidth bounding the middle 50% of the 
cumulative integrated power. 



To understand the physical mechanism controlling 
the steady and unsteady states we examine the power 
spectra of D{t) in the glycerol solution experiments. 
Figures [TWa) and (b) represent the power spectra for 
the steady (+) and unsteady (•) perturbation behavior 
shown in Figure [TTT a). In the steady case the funda- 
mental frequency (/o), which is the first harmonic of the 
power spectra, represents the rate at which perturbations 
form along the fiber (e.g., fo — 14.45 perturbations/s in 
the experiment shown in Figure E^b)); in the unsteady 
case the fundamental peak is much broader so that fo is 
less well defined. As a function of increasing flow rate, the 
fundamental frequency (i.e., rate of perturbation forma- 
tion) increases linearly when the perturbation dynamics 
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FIG. 14: Interquartile Region (IQR), or frequency bandwidth, 
as a function of flow rate for the glycerol solution experiments. 
An abrupt transition in the bandwidth occurs at the transi- 
tion flow rate, Q c - 



are steady (Q < Q c ), and is scattered about w 15 per- 
turbations/s when the dynamics are unsteady (Q > Q c ) 
(see Figure [T37a)). Another feature in the power spec- 
tra distinguishing steady and unsteady behavior is in the 
frequency bandwidth supporting the fundamental peak; 
the bandwidth of the unsteady spectra is larger than 
the steady spectra in Figure [T3J We characterize the 
bandwidth of the fundamental peak by measuring the 
Interquartile Region (IQR). The IQR is defined as the 
frequency bandwidth bounding the middle 50% of the 
(normalized) cumulative integrated power under the fun- 
damental peak; an example is shown in Figure [P37 b) cor- 
responding to the power spectra in Figure [F2T b) [3 II]. The 
IQR, or bandwidth, measures the modulation of the fun- 
damental frequency, or more physically, the modulation 
of the rate at which perturbations form along the fiber. A 
jump in the bandwidth occurs at the transition flow rate 
Q c in the glycerol solution experiments, as shown in Fig- 
ure [Mj For Q < Q c , the bandwidth is nearly zero, thus 
the rate of perturbation formation is nearly constant re- 
sulting in longer-time steady perturbation behavior. For 
Q > Qc the bandwidth is sizable and increases with 
increasing flow rate, thus there is a significant modula- 
tion of the rate at which perturbations form. It is this 
large modulation that results in the longer-time unsteady 
dynamics of the perturbations. While the transition in 
Figure [14] is striking, it is not entirely clear whether it is 
a subcritical or supercritical transition, and if subcritical, 
whether the transition is hysteretic. 



V. CONCLUSIONS 

In an experimental study, we examine the motion of 
an annular viscous film flowing under the influence of 
gravity down the outside of a vertical fiber. We find the 
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unperturbed flow is well approximated by a steady, uni- 
directional parallel flow when Re < 10. The dynamics 
of the perturbed flow can be divided into three stages: 
(i) initial exponential growth of the perturbation am- 
plitude accompanied by a decrease in the perturbation 
wavelength; (ii) nonlinear saturation of the perturbation 
amplitude and wavelength; and (iii) longer-time behavior 
in which the perturbation wavelength may (unsteady) or 
may not (steady) vary along the film. During the first 
stage, we find linear stability theory results developed 
from a long-wave Stokes flow model [15j are in excellent 
agreement with the initial growth of perturbations mea- 
sured in experiments. The agreement between linear sta- 
bility results developed from a moderate Reynolds num- 
ber model [IH and experimental data are not as strong 
as in the Stokes flow case. 

A close examination of the longer-time steady and un- 
steady behavior of intcrfacial perturbations is shown to 
be correlated to the range of: (i) the rate of exponential 
growth of the perturbation amplitude; and (ii) the lo- 
cation along the fiber where perturbations initially form. 
In particular, we find the rate of growth of the amplitude 
and the location along the fiber where perturbations form 
is nearly constant for the steady case, and varies over a 
range of values in the unsteady case. Furthermore, we 
find the transition in the longer-time perturbation dy- 
namics from unsteady to steady behavior at a critical 
flow rate occurs because of a transition in the rate at 
which perturbations naturally form along the free surface 



of the film. In the steady case, the rate of perturbation 
formation is nearly constant, resulting in the perturba- 
tions remaining equally spaced as they travel with the 
same terminal speed down the fiber. In the unsteady 
case, the rate of perturbation formation is modulated 
which results in the modulation of the initial speed and 
spacing between perturbations and ultimately leads to 
the coalescence of perturbations further down the fiber. 
It is not clear whether this transition is subcritical or su- 
percritical, and if subcritical, whether the transition is 
hysteretic. 
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APPENDIX 

Formula for the constant coefficients described in the 
dispersion relation derived by Sisoev et al. (l8j are: 



where 



«o,o = (l + §) (3 + fr), ao,3 = ^ (l + |) , (A.la) 

^ = I^w(^ 201 - 6 )-i( 1 + ^)(^T7)) 2 ' (A ' lb) 

ai ' x = iwmr ai,0 = i( 1+ i) ! (A - lc) 



b = 2(l + e) 6 (ln(l + e)) 2 + (2-3(l + e) 2 )(l+e) 4 ln(l + e) 

+ ^(l + e) 6 -^l + e) 4 + ^(l + £ ) 2 -i, (A.ld) 

h = l^(l + £ f -10 £ (l + e ) 3 + |(l + £ )+4 £ (l + £ ) 5 (ln(l + e)) 2 (3 +i ^^ 

+ e(2 - 3(1 + em + ef ln(l + s) (£- - ^|±^ + — ^ ^ — ) , (A.le) 
e(l + e) 3 (41n(l + e) + 1) - e - 7e 2 - 9e 3 - 3 £ 4 

1,2.^3 3,4 3 ' ( A - lf ) 



(1 + e) 4 ln(l + e)-e- \e 2 - 3e 3 
6(f) = 4(1 + £ ) 4 ln(l + e) ((3(1 + £ ) 2 - !)((! + ef 1)) 
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